Priming avocado with sodium hydrosulfide prior to frost conditions induces the expression of genes involved in protection and stress responses.

Priming plants with chemical agents has been extensively investigated as a means for improving their tolerance to many biotic and abiotic stresses. Earlier, we showed that priming young avocado (Persea americana Mill cv. 'Hass') trees with sodium hydrosulfide (NaHS), a donor of hydrogen sulfide, improves the response of photosynthesis to simulated frost (cold followed by high light) conditions. In the current study, we performed a transcriptome analysis to gain insight into the molecular response of avocado 'Hass' leaves to frost, with or without NaHS priming. The analysis revealed 2144 (down-regulated) and 2064 (up-regulated) differentially expressed genes (DEGs) common to both non-primed and primed trees. Non-primed trees had 697 (down) and 559 (up) unique DEGs, while primed trees exhibited 1395 (down) and 1385 (up) unique DEGs. We focus on changes in the expression patterns of genes encoding proteins involved in photosynthesis, carbon cycle, protective functions, biosynthesis of isoprenoids and abscisic acid (ABA), as well as ABA-regulated genes. Notably, the differential expression results depict the enhanced response of primed trees to the frost and highlight gene expression changes unique to primed trees. Amongst these are up-regulated genes encoding pathogenesis-related proteins, heat shock proteins, enzymes for ABA metabolism, and ABA-induced transcription factors. Extending the priming experiments to field conditions, which showed a benefit to the physiology of trees following chilling, suggests that it can be a possible means to improve trees' response to cold stress under natural winter conditions.

Phenylalanine treatment induces tomato resistance to Tuta absoluta via increased accumulation of benzenoid/phenylpropanoid volatiles serving as defense signals.

Tuta absoluta ("leafminer"), is a major pest of tomato crops worldwide. Controlling this insect is difficult due to its efficient infestation, rapid proliferation, and resilience to changing weather conditions. Furthermore, chemical pesticides have only a short-term effect due to rapid development of T. absoluta strains. Here, we show that a variety of tomato cultivars, treated with external phenylalanine solutions exhibit high resistance to T. absoluta, under both greenhouse and open field conditions, at different locations. A large-scale metabolomic study revealed that tomato leaves absorb and metabolize externally given Phe efficiently, resulting in a change in their volatile profile, and repellence of T. absoluta moths. The change in the volatile profile is due to an increase in three phenylalanine-derived benzenoid phenylpropanoid volatiles (BPVs), benzaldehyde, phenylacetaldehyde, and 2-phenylethanol. This treatment had no effect on terpenes and green leaf volatiles, known to contribute to the fight against insects. Phe-treated plants also increased the resistance of neighboring non-treated plants. RNAseq analysis of the neighboring non-treated plants revealed an exclusive upregulation of genes, with enrichment of genes related to the plant immune response system. Exposure of tomato plants to either benzaldehyde, phenylacetaldehyde, or 2-phenylethanol, resulted in induction of genes related to the plant immune system that were also induced due to neighboring Phe-treated plants. We suggest a novel role of phenylalanine-derived BPVs as mediators of plant-insect interactions, acting as inducers of the plant defense mechanisms.

Deprivation of Sexual Reproduction during Garlic Domestication and Crop Evolution.

Garlic, originating in the mountains of Central Asia, has undergone domestication and subsequent widespread introduction to diverse regions. Human selection for adaptation to various climates has resulted in the development of numerous garlic varieties, each characterized by specific morphological and physiological traits. However, this process has led to a loss of fertility and seed production in garlic crops. In this study, we conducted morpho-physiological and transcriptome analyses, along with whole-genome resequencing of 41 garlic accessions from different regions, in order to assess the variations in reproductive traits among garlic populations. Our findings indicate that the evolution of garlic crops was associated with mutations in genes related to vernalization and the circadian clock. The decline in sexual reproduction is not solely attributed to a few mutations in specific genes, but is correlated with extensive alterations in the genetic regulation of the annual cycle, stress adaptations, and environmental requirements. The regulation of flowering ability, stress response, and metabolism occurs at both the genetic and transcriptional levels. We conclude that the migration and evolution of garlic crops involve substantial and diverse changes across the entire genome landscape. The construction of a garlic pan-genome, encompassing genetic diversity from various garlic populations, will provide further insights for research into and the improvement of garlic crops.

Promoter insertion leads to polyembryony in mango - a case of convergent evolution with citrus.

Sexual reproduction in plants is the main pathway for creating new genetic combinations in modern agriculture. In heterozygous plants, after the identification of a plant with desired traits, vegetative propagation (cloning) is the primary path to create genetically uniform plants. Another natural plant mechanism that creates genetically uniform plants (clones) is apomixis. In fruit crops like citrus and mango, sporophytic apomixis results in polyembryony, where seeds contain multiple embryos, one of which is sexually originated and the others are vegetative clones of the parent mother tree. Utilizing the mango genome and genetic analysis of a diverse germplasm collection, we identified MiRWP as the gene that causes polyembryony in mango. There is a strong correlation between a specific insertion in the gene's promoter region and altered expression in flowers and developing fruitlets, inducing multiple embryos. The MiRWP gene is an ortholog of CitRWP that causes polyembryony in citrus. Based on the data, we speculate that promoter insertion events, which occurred independently in citrus and mango, induced nucellar embryogenesis. The results suggest convergent evolution of polyembryony in the two species. Further work is required to demonstrate the utility of these genes (mango and citrus) in other biological systems as a tool for the clonal production of other crops.

Pathogen Eradication in Garlic in the Phytobiome Context: Should We Aim for Complete Cleaning?

Global food production is challenged by plant pathogens that cause significant crop losses. Fungi, bacteria, and viruses have long threatened sustainable and profitable agriculture. The danger is even higher in vegetatively propagated horticultural crops, such as garlic. Currently, quarantine, rouging infected plants, and control of natural vectors are used as the main means of disease and pest control in garlic crops. Agricultural biotechnology, meristem-tip culture, and cryotherapy offer solutions for virus eradication and for the multiplication of 'clean stocks', but at the same time, impact the symbiotic and beneficial components of the garlic microbiome. Our research involves the first metatranscriptomic analysis of the microbiome of garlic bulb tissue, PCR analyses, and a biological assay of endophytes and pathogens. We have demonstrated that in vitro sanitation methods, such as shoot tip culture or cryotherapy can alter the garlic microbiome. Shoot tip culture proved ineffective in virus elimination, but reduced bacterial load and eliminated fungal infections. Conversely, cryotherapy was efficient in virus eradication but demolished other components of the garlic microbiome. Garlic plants sanitized by cryotherapy exhibited a lower survival rate, and a longer in vitro regeneration period. The question arises whether total eradication of viruses, at the expense of other microflora, is necessary, or if a partial reduction in the pathogenic load would suffice for sanitized garlic production. We explore this question from both scientific and commercial perspectives.

Carbohydrate dynamics in Populus trees under drought: An expression atlas of genes related to sensing, translocation, and metabolism across organs.

In trees, nonstructural carbohydrates (NSCs) serve as long-term carbon storage and long-distance carbon transport from source to sink. NSC management in response to drought stress is key to our understanding of drought acclimation. However, the molecular mechanisms underlying these processes remain unclear. By combining a transcriptomic approach with NSC quantification in the leaves, stems, and roots of Populus alba under drought stress, we analyzed genes from 29 gene families related to NSC signaling, translocation, and metabolism. We found starch depletion across organs and accumulation of soluble sugars (SS) in the leaves. Activation of the trehalose-6-phosphate/SNF1-related protein kinase (SnRK1) signaling pathway across organs via the suppression of class I TREHALOSE-PHOSPHATE SYNTHASE (TPS) and the expression of class II TPS genes suggested an active response to drought. The expression of SnRK1α and ß subunits, and SUCROSE SYNTHASE6 supported SS accumulation in leaves. The upregulation of active transporters and the downregulation of most passive transporters implied a shift toward active sugar transport and enhanced regulation over partitioning. SS accumulation in vacuoles supports osmoregulation in leaves. The increased expression of sucrose synthesis genes and reduced expression of sucrose degradation genes in the roots did not coincide with sucrose levels, implying local sucrose production for energy. Moreover, the downregulation of invertases in the roots suggests limited sucrose allocation from the aboveground organs. This study provides an expression atlas of NSC-related genes that respond to drought in poplar trees, and can be tested in tree improvement programs for adaptation to drought conditions.

Aneuploidy Formation in the Filamentous Fungus Aspergillus flavus in Response to Azole Stress.

Aspergillus flavus is a mycotoxigenic fungus that contaminates many important agricultural crops with aflatoxin B1, the most toxic and carcinogenic natural compound. This fungus is also the second leading cause of human invasive aspergillosis, after Aspergillus fumigatus, a disease that is particularly prevalent in immunocompromised individuals. Azole drugs are considered the most effective compounds in controlling Aspergillus infections both in clinical and agricultural settings. Emergence of azole resistance in Aspergillus spp. is typically associated with point mutations in cyp51 orthologs that encode lanosterol 14α-demethylase, a component of the ergosterol biosynthesis pathway that is also the target of azoles. We hypothesized that alternative molecular mechanisms are also responsible for acquisition of azole resistance in filamentous fungi. We found that an aflatoxin-producing A. flavus strain adapted to voriconazole exposure at levels above the MIC through whole or segmental aneuploidy of specific chromosomes. We confirm a complete duplication of chromosome 8 in two sequentially isolated clones and a segmental duplication of chromosome 3 in another clone, emphasizing the potential diversity of aneuploidy-mediated resistance mechanisms. The plasticity of aneuploidy-mediated resistance was evidenced by the ability of voriconazole-resistant clones to revert to their original level of azole susceptibility following repeated transfers on drug-free media. This study provides new insights into mechanisms of azole resistance in a filamentous fungus. IMPORTANCE Fungal pathogens cause human disease and threaten global food security by contaminating crops with toxins (mycotoxins). Aspergillus flavus is an opportunistic mycotoxigenic fungus that causes invasive and noninvasive aspergillosis, diseases with high rates of mortality in immunocompromised individuals. Additionally, this fungus contaminates most major crops with the notorious carcinogen, aflatoxin. Voriconazole is the drug of choice to treat infections caused by Aspergillus spp. Although azole resistance mechanisms have been well characterized in clinical isolates of Aspergillus fumigatus, the molecular basis of azole resistance in A. flavus remains unclear. Whole-genome sequencing of eight voriconazole-resistant isolates revealed that, among other factors, A. flavus adapts to high concentrations of voriconazole by duplication of specific chromosomes (i.e., aneuploidy). Our discovery of aneuploidy-mediated resistance in a filamentous fungus represents a paradigm shift, as this type of resistance was previously thought to occur only in yeasts. This observation provides the first experimental evidence of aneuploidy-mediated azole resistance in the filamentous fungus A. flavus.

Non-Aqueous Isolation and Enrichment of Glandular Capitate Stalked and Sessile Trichomes from Cannabis sativa.

This paper presents a protocol for the convenient and high-throughput isolation and enrichment of glandular capitate stalked and sessile trichomes from Cannabis sativa. The biosynthetic pathways for cannabinoid and volatile terpene metabolism are localized primarily in the Cannabis trichomes, and isolated trichomes are beneficial for transcriptome analysis. The existing protocols for isolating glandular trichomes for transcriptomic characterization are inconvenient and deliver compromised trichome heads and a relatively low amount of isolated trichomes. Furthermore, they rely on expensive apparatus and isolation media containing protein inhibitors to avoid RNA degradation. The present protocol suggests combining three individual modifications to obtain a large amount of isolated glandular capitate stalked and sessile trichomes from C. sativa mature female inflorescences and fan leaves, respectively. The first modification involves substituting liquid nitrogen for the conventional isolation medium to facilitate the passage of trichomes through the micro-sieves. The second modification involves using dry ice to detach the trichomes from the plant source. The third modification involves passing the plant material consecutively through five micro-sieves of diminishing pore sizes. Microscopic imaging demonstrated the effectiveness of the isolation technique for both trichome types. In addition, the quality of RNA extracted from the isolated trichomes was appropriate for downstream transcriptomic analysis.

Discovery of an unrecognized nidovirus associated with granulomatous hepatitis in rainbow trout.

Rainbow trout (Oncorhynchus mykiss) is the principal species of inland-farmed fish in the Western hemisphere. Recently, we diagnosed in farmed rainbow trout a disease in which the hallmark is granulomatous-like hepatitis. No biotic agents could be isolated from lesions. Still, unbiased high-throughput sequencing and bioinformatics analyses revealed the presence of a novel piscine nidovirus that we named "Trout Granulomatous Virus" (TGV). TGV genome (28,767 nucleotides long) is predicted to encode non-structural (1a and 1 ab) and structural (S, M, and N) proteins that resemble proteins of other known piscine nidoviruses. High loads of TGV transcripts were detected by quantitative RT-PCR in diseased fish and visualized in hepatic granulomatous sites by fluorescence in situ hybridization. Transmission electron microscopy (TEM) revealed coronavirus-like particles in these lesions. Together, these analyses corroborated the association of TGV with the lesions. The identification and detection of TGV provide means to control TGV spread in trout populations.

Lysine 27 of histone H3.3 is a fine modulator of developmental gene expression and stands as an epigenetic checkpoint for lignin biosynthesis in Arabidopsis.

Chromatin is a dynamic platform within which gene expression is controlled by epigenetic modifications, notably targeting amino acid residues of histone H3. Among them is lysine 27 of H3 (H3K27), the trimethylation of which by the Polycomb Repressive Complex 2 (PRC2) is instrumental in regulating spatiotemporal patterns of key developmental genes. H3K27 is also subjected to acetylation and is found at sites of active transcription. Most information on the function of histone residues and their associated modifications in plants was obtained from studies of loss-of-function mutants for the complexes that modify them. To decrypt the genuine function of H3K27, we expressed a non-modifiable variant of H3 at residue K27 (H3.3K27A ) in Arabidopsis, and developed a multi-scale approach combining in-depth phenotypical and cytological analyses, with transcriptomics and metabolomics. We uncovered that the H3.3K27A variant causes severe developmental defects, part of them are reminiscent of PRC2 mutants, part of them are new. They include early flowering, increased callus formation and short stems with thicker xylem cell layer. This latest phenotype correlates with mis-regulation of phenylpropanoid biosynthesis. Overall, our results reveal novel roles of H3K27 in plant cell fates and metabolic pathways, and highlight an epigenetic control point for elongation and lignin composition of the stem.

Two independent loss-of-function mutations in anthocyanidin synthase homeologous genes are responsible for the all-green phenotype of sweet basil.

Sweet basil, Ocimum basilicum L., is an important culinary herb grown worldwide. Although basil is green, many landraces, breeding lines, and exotic cultivars have purple stems and flowers. This anthocyanin pigmentation is unacceptable in traditional Italian basil used for Pesto sauce production. In the current study, we aimed to resolve the genetics that underlines the different colors. We used the recently published sweet basil genome to map quantitative trait loci (QTL) for flower and stem color in a bi-parental F2 population. It was found that the pigmentation is governed by a single QTL, harboring an anthocyanidin synthase (ANS) gene (EC 1.14.20.4). Further analysis revealed that the basil genome harbors two homeologous ANS genes, each carrying a loss-of-function mutation. ObANS1 carries a single base pair insertion resulting in a frameshift, and ObANS2 carries a missense mutation within the active site. In the purple-flower parent, ANS1 is functional, and ANS2 carries a nonsense mutation. The functionality of the ObANS1 active allele was validated by complementation assay in an Arabidopsis ANS mutant. Moreover, we have restored the functionality of the missense-mutated ObANS2 using site-directed activation. We found that the non-functional alleles were expressed to similar levels as the functional allele, suggesting polyploids invest futile effort in expressing non-functional genes, offsetting their advantageous redundancy. This work demonstrated the usefulness of the genomics and genetics of basil to understand the basic mechanism of metabolic traits and raise fundamental questions in polyploid plant biology.

A framework for the targeted recruitment of crop-beneficial soil taxa based on network analysis of metagenomics data.

BACKGROUND: The design of ecologically sustainable and plant-beneficial soil systems is a key goal in actively manipulating root-associated microbiomes. Community engineering efforts commonly seek to harness the potential of the indigenous microbiome through substrate-mediated recruitment of beneficial members. In most sustainable practices, microbial recruitment mechanisms rely on the application of complex organic mixtures where the resources/metabolites that act as direct stimulants of beneficial groups are not characterized. Outcomes of such indirect amendments are unpredictable regarding engineering the microbiome and achieving a plant-beneficial environment. RESULTS: This study applied network analysis of metagenomics data to explore amendment-derived transformations in the soil microbiome, which lead to the suppression of pathogens affecting apple root systems. Shotgun metagenomic analysis was conducted with data from 'sick' vs 'healthy/recovered' rhizosphere soil microbiomes. The data was then converted into community-level metabolic networks. Simulations examined the functional contribution of treatment-associated taxonomic groups and linked them with specific amendment-induced metabolites. This analysis enabled the selection of specific metabolites that were predicted to amplify or diminish the abundance of targeted microbes functional in the healthy soil system. Many of these predictions were corroborated by experimental evidence from the literature. The potential of two of these metabolites (dopamine and vitamin B12) to either stimulate or suppress targeted microbial groups was evaluated in a follow-up set of soil microcosm experiments. The results corroborated the stimulant's potential (but not the suppressor) to act as a modulator of plant beneficial bacteria, paving the way for future development of knowledge-based (rather than trial and error) metabolic-defined amendments. Our pipeline for generating predictions for the selective targeting of microbial groups based on processing assembled and annotated metagenomics data is available at https://github.com/ot483/NetCom2 . CONCLUSIONS: This research demonstrates how genomic-based algorithms can be used to formulate testable hypotheses for strategically engineering the rhizosphere microbiome by identifying specific compounds, which may act as selective modulators of microbial communities. Applying this framework to reduce unpredictable elements in amendment-based solutions promotes the development of ecologically-sound methods for re-establishing a functional microbiome in agro and other ecosystems. Video Abstract.

An alternative pathway to plant cold tolerance in the absence of vacuolar invertase activity.

To cope with cold stress, plants have developed antioxidation strategies combined with osmoprotection by sugars. In potato (Solanum tuberosum) tubers, which are swollen stems, exposure to cold stress induces starch degradation and sucrose synthesis. Vacuolar acid invertase (VInv) activity is a significant part of the cold-induced sweetening (CIS) response, by rapidly cleaving sucrose into hexoses and increasing osmoprotection. To discover alternative plant tissue pathways for coping with cold stress, we produced VInv-knockout lines in two cultivars. Genome editing of VInv in 'Désirée' and 'Brooke' was done using stable and transient expression of CRISPR/Cas9 components, respectively. After storage at 4°C, sugar analysis indicated that the knockout lines showed low levels of CIS and maintained low acid invertase activity in storage. Surprisingly, the tuber parenchyma of vinv lines exhibited significantly reduced lipid peroxidation and reduced H2 O2 levels. Furthermore, whole plants of vinv lines exposed to cold stress without irrigation showed normal vigor, in contrast to WT plants, which wilted. Transcriptome analysis of vinv lines revealed upregulation of an osmoprotectant pathway and ethylene-related genes during cold temperature exposure. Accordingly, higher expression of antioxidant-related genes was detected after exposure to short and long cold storage. Sugar measurements showed an elevation of an alternative pathway in the absence of VInv activity, raising the raffinose pathway with increasing levels of myo-inositol content as a cold tolerance response.

Comparing adventitious root-formation and graft-unification abilities in clones of Argania spinosa.

Argania spinosa trees have attracted attention in recent years due to their high resistance to extreme climate conditions. Initial domestication activities practiced in Morocco. Here we report on selection and vegetative propagation of A. spinosa trees grown in Israel. Trees yielding relatively high amounts of fruit were propagated by rooting of stem cuttings. High variability in rooting ability was found among the 30 clones selected. In-depth comparison of a difficult-to-root (ARS7) and easy-to-root (ARS1) clone revealed that the rooted cuttings of ARS7 have a lower survival rate than those of ARS1. In addition, histological analysis of the adventitious root primordia showed many abnormal fused primordia in ARS7. Hormone profiling revealed that while ARS1 accumulates more cytokinin, ARS7 accumulates more auxin, suggesting different auxin-to-cytokinin ratios underlying the different rooting capabilities. The hypothesized relationship between rooting and grafting abilities was addressed. Reciprocal grafting was performed with ARS1/ARS7 but no significant differences in the success of graft unification between the trees was detected. Accordingly, comparative RNA sequencing of the rooting and grafting zones showed more differentially expressed genes related to rooting than to grafting between the two trees. Clustering, KEGG and Venn analyses confirmed enrichment of genes related to auxin metabolism, transport and signaling, cytokinin metabolism and signaling, cell wall modification and cell division in both regions. In addition, the differential expression of some key genes in ARS1 vs. ARS7 rooting zones was revealed. Taken together, while both adventitious root-formation and graft-unification processes share response to wounding, cell reprogramming, cell division, cell differentiation and reconnection of the vasculature, there are similar, but also many different genes regulating the two processes. Therefore an individual genotype can have low rooting capacity but good graft-unification ability.

Integrated Genomic and Transcriptomic Elucidation of Flowering in Garlic.

Commercial cultivars of garlic are sterile, and therefore efficient breeding of this crop is impossible. Recent restoration of garlic fertility has opened new options for seed production and hybridization. Transcriptome catalogs were employed as a basis for garlic genetic studies, and in 2020 the huge genome of garlic was fully sequenced. We provide conjoint genomic and transcriptome analysis of the regulatory network in flowering garlic genotypes. The genome analysis revealed phosphatidylethanolamine-binding proteins (PEBP) and LEAFY (LFY) genes that were not found at the transcriptome level. Functions of TFL-like genes were reduced and replaced by FT-like homologs, whereas homologs of MFT-like genes were not found. The discovery of three sequences of LFY-like genes in the garlic genome and confirmation of their alternative splicing suggest their role in garlic florogenesis. It is not yet clear whether AsLFY1 acts alone as the "pioneer transcription factor" or AsLFY2 also provides these functions. The presence of several orthologs of flowering genes that differ in their expression and co-expression network advocates ongoing evolution in the garlic genome and diversification of gene functions. We propose that the process of fertility deprivation in garlic cultivars is based on the loss of transcriptional functions of the specific genes.

Non-photoperiodic transition of female cannabis seedlings from juvenile to adult reproductive stage.

KEY MESSAGE: Vegetative-to-reproductive phase transition in female cannabis seedlings occurs autonomously with the de novo development of single flowers. To ensure successful sexual reproduction, many plant species originating from seedlings undergo juvenile-to-adult transition. This phase transition precedes and enables the vegetative-to-reproductive shift in plants, upon perception of internal and/or external signals such as temperature, photoperiod, metabolite levels, and phytohormones. This study demonstrates that the juvenile seedlings of cannabis gradually shift to the adult vegetative stage, as confirmed by the formation of lobed leaves, and upregulation of the phase-transition genes. In the tested cultivar, the switch to the reproductive stage occurs with the development of a pair of single flowers in the 7th node. Histological analysis indicated that transition to the reproductive stage is accomplished by the de novo establishment of new flower meristems which are not present in a vegetative stage, or as dormant meristems at nodes 4 and 6. Moreover, there were dramatic changes in the transcriptomic profile of flowering-related genes among nodes 4, 6, and 7. Downregulation of flowering repressors and an intense increase in the transcription of phase transition-related genes occur in parallel with an increase in the transcription of flowering integrators and meristem identity genes. These results support and provide molecular evidence for previous findings that cannabis possesses an autonomous flowering mechanism and the transition to reproductive phase is controlled in this plant mainly by internal signals.

Ripening of Pomegranate Skin as Revealed by Developmental Transcriptomics.

The appearance of pomegranate (Punica granatum L.) fruit is highly important for its marketing. The primary concerns are obtaining sufficient red pigment accumulation and minimal cracking of the fruit skin (the outer red layer of the peel). We analyzed the skin transcriptome of pomegranate cv. Wonderful at distinct time points of fruit development to characterize the processes that occur in the skin during fruit ripening and which may reflect on processes in the whole fruit, such as the non-climacteric nature of pomegranate. The data suggested a ripening mechanism in pomegranate skin that differs from that in strawberry-the model plant for non-climacteric fruit where abscisic acid is the growth regulator that drives ripening-involving ethylene, polyamine, and jasmonic acid pathways. The biosynthetic pathways of important metabolites in pomegranate-hydrolyzable tannins and anthocyanins-were co-upregulated at the ripening stage, in line with the visual enhancement of red coloration. Interestingly, cuticle- and cell-wall-related genes that showed differential expression between the developmental stages were mainly upregulated in the skin of early fruit, with lower expression at mid-growth and ripening stages. Nevertheless, lignification may be involved in skin hardening in the mature fruit.

The SWC4 subunit of the SWR1 chromatin remodeling complex is involved in varying virulence of Metarhizium brunneum isolates offering role of epigenetic regulation of pathogenicity.

The host - pathogen interaction is a multifactorial process subject to a co-evolutionary arms race consisting of rapid changes in both host and pathogen, controlled at the genetic and epigenetic levels. Previously, we showed intra-species variation in disease progression and pathogenicity in aphids for Metarhizium brunneum isolates MbK and Mb7. Herein, we compared genomic, epigenetic, and metabolomic variations between these isolates and their effects on pathogenicity. Genomic variation could not completely explain the observed differences between the isolates. However, differential N6-adenine methylation (6 mA) and its correlation to reduced expression of the essential SWC4 subunit of SWR1 chromatin-remodelling complex (SWR1-C) led us to hypothesize a role for swc4 in the varying pathogenicity. Mutagenesis of the essential swc4 gene in MbKisolate resulted in reduction of secondary-metabolite (SM) secretion and impaired virulence in Galleria mellonella. Our results suggest the role of SWC4 in the regulation of SMs and the role of both SWC4 and SWR1-C in virulence of M. brunneum isolates. A better understanding of epigenetic regulation of SM production and secretion in entomopathogenic fungi may enable theirmanipulation for better biocontrol performance, and expand possibilities for environmentally friendly pest control.